Precise measurement of the CP violation parameter sin20i in B° — > (cc)K° decays 
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We present a precise measurement of the CP violation parameter sin 20i and the direct CP 
violation parameter Af using the final data sample of 772 x fO B BB pairs collected at the T(45) 
resonance with the Belle detector at the KEKB asymmetric-energy e + e~ collider. One neutral 
B meson is reconstructed in a J/ipKl, ip{2S)K%, XciK° s or J/ipK° L CP-eigenstate and its flavor 
is identified from the decay products of the accompanying B meson. From the distribution of 
proper time intervals between the two B decays, we obtain the following CP violation parameters: 
sin 20i = 0.667 ± 0.023(stat) ± 0.012(syst) and Af = 0.006 ± 0.016(stat) ± 0.012(syst). 

PACS numbers: 11.30.Er, 12.15.Hh, 13.25.Hw 



In the standard model (SM), CP violation in the quark 
sector is described by the Kobayashi-Maskawa (KM) the- 
ory [l[ in which the quark-mixing matrix has a single irre- 
ducible complex phase that gives rise to all CP-violating 
asymmetries. In the decay chain T(AS) — > B°B° — > 
/cp/tag, where one of the B mesons decays at time tcp 
to a CP-eigenstate fc p and the other B meson decays at 
time itag to a final state /t ag that distinguishes between 
B° and B°, the decay rate has a time dependence in the 
T(4S) rest frame 0] given by 



V{At) 



,-\At\/T B o 

4r B o 



Sf sm(Am d At) 



A f cos(Am d At) \. (1) 



Here Sf and Af are CP violation parameters, t b o is 
the B° lifetime, Am d is the mass difference between the 
two neutral B mass eigenstates, At = tcp — itag, and 
the b- flavor charge q = +1 (—1) when the tagging B 
meson is a B° (B°). With very small theoretical uncer- 
tainty 0, the SM predicts Sf = — £/ sin20i and Af = 
for the b —> ccs transition, where £/ = +1 (—1) cor- 
responds to CP-even (-odd) final states and 4>i is an 
interior angle of the KM unitarity triangle, defined as 
01 = ^g[-V cd V* b /V td V* b } 1. The BaBar and Belle 



collaborations have published several determinations of 
sin 20i since the first observation [3, [B[ ; previous results 
used 465 x 10 6 @ and 535 x 10 6 BB pairs, respectively 

With recently available experimental results, not only 
sin 20i but also other measurements of the sides of the 
unitarity triangle and other CP violation measurements 
make it possible to test the consistency of the KM 
scheme. The indirect determination of the angle 0i de- 
viates by 2.7cr from the current world average for the 
direct determination of sin 20i [8( . Equivalently, the 
B^ — > r^Vr branching fraction and the resulting value 
of \Vub\ differ by 2.8ct from the prediction of the global 
fit Q, where the sin20i value gives the most stringent 
constraint on the indirect measurement. Furthermore, 
time-dependent CP violation in the neutral B meson de- 
cays mediated by flavor-changing b — > s transitions may 
deviate from CP violation in the b — > ccs case because of 
possible additional quantum loops [9] . To clarify whether 
new physics contributes to CP-violating phenomena or 
B — > t^v t decays, it is very important to determine 
sin 20i, the SM reference, as precisely as possible. 

In this Letter, we describe the final Belle measure- 
ment of sin 20i and Af in b —> ccs induced B decays to 
fcp- The B decays to the CP-odd eigenstates, fcp = 
J/yjK° s , vb(2S)K% and XciK%, and the CP-even eigen- 
state, fcp = J li>Kj\, are reconstructed using 772 x 10 6 
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BB pairs, the entire data sample accumulated on the 
Y(4S) resonance with the Belle detecto r flOi] at the 
KEKB asymmetric-energy e + e~ collider 11]. Two in- 
ner detector configurations were used. A 2.0 cm radius 
beampipe and a 3-layer silicon vertex detector (SVD) 
were used for the first data sample that contains 152 x 10 6 
BB pairs. The remaining 620 x 10 6 BB pairs were accu- 
mulated with a 1.5 cm radius beampipe, a 4-layer silicon 
vertex detector and a small-cell inner drift chamber. The 
latter data sample has been recently reprocessed using a 
new charged track reconstruction algorithm, which sig- 
nificantly increased the reconstruction efficiency for the 
B° — > (cc)Kg decay modes. In particular, the gain for 
the B° -> J/i}jK% decay mode is 18%. 

The T(45) is produced with a Lorentz boost of (3j — 
0.425 nearly along the z-axis, which is antiparallel to the 
positron beam direction. Since the B° and B° mesons are 
approximately at rest in the T(4S I ) center-of-mass system 
(CM), At can be determined from the displacement in z 
between the fcp and /t ag decay vertices: At ~ (zcp — 
z tag )/(/3 7 c) = Az/(/3 7 c). 

Charged tracks reconstructed in the central drift cham- 
ber (CDC), except for tracks from Kg —> ir + ir~ decays, 
are required to originate from the interaction point (IP). 
We distinguish charged kaons from pions based on a kaon 
(pion) likelihood Cjcm derived from the time-of-flight 
scintillation counters, aerogel threshold Cherenkov coun- 
ters (ACC), and dE/dx measurements in the CDC. Elec- 
tron identification is based on the ratio of the electro- 
magnetic calorimeter (ECL) cluster energy to the parti- 
cle momentum as well as a combination of dE/dx mea- 
surements in the CDC, the ACC response, and the posi- 
tion and shape of the electromagnetic shower. Muons are 
identified by track penetration depth and hit scatter in 
the muon detector (KLM). Photons are identified as iso- 
lated ECL clusters that are not matched to any charged 
track. 

For the J/ipK° s , J/ipKl and ip(2S)K s modes, event 
selection is the same as in our previous analyses 0, [H| , 
where J/tp mesons are reconstructed via their decays 
to £+£- (£ = e,/x) and the ip(2S) mesons to £ + £" or 
J/ipir + ir~. For the modes J/i/jK^ and XciK s , in which 
the Xci is reconstructed in the J/ipj final state, both 
J I i/j daughter tracks must be positively identified as lep- 
tons, whereas for the J/ipK s and ip(2S)K s modes, at 
least one daughter must satisfy this requirement. Any 
other track having an ECL energy deposit consistent 
with a minimum ionizing particle is accepted as a muon 
candidate and any track satisfying either the dE/dx 
or the ECL shower energy requirements is retained as 
an electron candidate. For J/ip — > e + e~ decays, the 
e charmonium daughters are combined with photons 
found within 50 mrad of the e + or e~ direction in or- 
der to account partially for final-state radiation and 
bremsstrahlung. In order to accommodate the remain- 
ing radiative tails, an asymmetric invariant mass require- 
ment is used to select J/tp and ip(2S) decays in dilep- 
ton modes, -150 MeV/c 2 < M e+e - - < 36 MeV/c 2 



and -60 MeV/c 2 < ^- - < 36 MeV/c 2 , where 
M^/j denotes either the nominal J/ip or ip(2S) mass. 
For ip(2S) — > J/ipTr + Tr~ candidates, we require a mass 
difference of 580 MeV/c 2 < M e + e - %+7r - - M t+l - < 
600 MeV/c 2 , and xa ~ > J/^Pl candidates are required 
to have a mass difference of 385.0 MeV/c 2 < M e + e - 7 — 
M e + e - < 430.5 MeV/c 2 . For each charmonium candi- 
date, vertex-constrained and mass-constrained fits are 
applied to improve its momentum resolution. 

Candidate K s —> ir + ir~ decays are selected by require- 
ments on their invariant mass, flight length and consis- 
tency between the K s momentum direction and vertex 
position. Candidate K® mesons are selected from ECL 
and/or KLM hit patterns that are consistent with the 
presence of a shower induced by a K® meson. The cen- 
troid of the K\ candidate shower is required to be within 
a 45° cone centered on the K\ direction calculated from 
the two-body B decay kinematics and the momentum of 
the reconstructed J/?/ 1 meson. 

For B — >• fcp candidate reconstruction in modes other 
than J/ijjKl, B candidates are identified by two kine- 
matic variables: the energy difference AE = E* B — 
^beam ana - the beam-energy constrained mass Mb c = 

V^bcam) 2 ~ CPs) 2 . where ^bcam the CM beam en- 
ergy, and E* B (p* B ) is the CM energy (momentum) of the 
reconstructed B candidate. The B° —t J/ipK® candi- 
dates are identified by the value of p* B calculated using a 
two-body decay kinematic assumption. 

The ^-flavor of the accompanying B meson is identi- 
fied from inclusive properties of particles that are not 
associated with the reconstructed B° — > fcp decay [l3j |. 
The tagging information is represented by two parame- 
ters, the 6-flavor charge q and purity r. The parameter 
r is an event- by- event, MC-determined flavor-tagging di- 
lution factor that ranges from r — for no flavor dis- 
crimination to r — 1 for unambiguous flavor assignment. 
The data are sorted into seven intervals of r. For events 
with r > 0.1, the wrong tag fractions for six r inter- 
vals, wi(l = 1,6), and their differences between B° and 
B° decays, Awi, are determined from semileptonic and 
hadronic b — » c decays [HI, If r < 0.1, the wrong 
tag fraction is set to 0.5, and therefore the tagging infor- 
mation is not used. The total effective tagging efficiency, 
E(/i x (1 - 2u;/) 2 ), is determined to be 0.298 ± 0.004, 
where // is the fraction of events in the category I. 

The vertex position for the fcp decay is reconstructed 
using J/tjj or tp(2S) daughter tracks that have a minimum 
number of SVD hits, while the /tag vertex is determined 
from well-reconstructed tracks that arc not assigned to 
fcp \XM ■ A constraint on the IP profile in the plane per- 
pendicular to the z-axis is used with the selected tracks. 
With this procedure, we are able to determine a vertex 
even in the case where only one track has sufficient asso- 
ciated SVD hits. The fractions of the single track vertices 
f° r I CP and / tag are about 12% and 23%, respectively. 

For a single track vertex, the estimated error of the z 
coordinate, a z , is the indicator of the vertex fit qual- 
ity and is required to be less than 500 /im. On the 
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other hand, a vertex reconstructed using two or more 
tracks is characterized by a more robust goodness-of- 
flt indicator. In the previous analysis Q, the value of 
X 2 of the vertex calculated solely along the z direction 
was used. This is now replaced by h, the value of x 2 
in three-dimensional space calculated using the charged 
tracks without using the interaction-region profile's con- 
straint [l5j |. A detailed MC study indicates that h is a 
superior indicator of the vertex goodness-of-fit because it 
is less sensitive to the specific B decay mode; in particu- 
lar, h shows a smaller mode dependence for the vertices 
reconstructed from B — > J/tpX and B — > D^X decays, 
which are used as control samples to determine the ver- 
tex resolution parameters. In the multiple-track vertex 
case, h < 50 and a z < 200 /im are required. For candi- 
date events in which both B vertices are reconstructed, 
we retain only those events where the B vertices satisfy 
|At| < 70 ps for further analysis. 

For the candidate events in which both flavor tag- 
ging and vertex reconstruction succeed, the signal yield 
and purity for each mode are obtained from an un- 
binned maximum-likelihood fit to the two-dimensional 
AE — Mbc distribution for fcp modes with a Kg me- 
son, and to the p* B distribution for J/ipK®. The back- 
ground mainly comes from BB events in which one of 
the B meson decays into a final state containing a cor- 
rectly reconstructed J/ip, i.e., the B — > J/ipX process. 
In order to determine this background distribution, a 
B J/tpX MC sample corresponding to 100 times the 
integrated luminosity of data is used. An estimate of 
other combinatorial backgrounds is obtained from the 
M e +£- sideband. For CP-odd modes, the signal dis- 
tribution is modeled with a Gaussian function in Mb c 
and a double Gaussian function in AE. The fits to 
determine signal yields for these modes are performed 
in the region 5.2 GeV/c 2 < M bc < 5.3 GeV/c 2 and 
-0.1 GeV < AE < 0.2 GeV. The p* B signal shape for 
J/ipK^ is determined from MC events. The requirement 
p* B < 2.0 GeV/c is used in the fit to estimate the signal 
yield as well as the contribution of three categories of 
background: those with a real (those that are correctly 
reconstructed) J/tp and a real K^, those with a real J/ip 
and a fake (those that are incorrectly reconstructed 
from electronic noise or electro- magnetic showers), and 
events with a fake J '/ip (those that are background com- 
binations). The Mbc distribution for a stringent AE re- 
quirement {\AE\ < 40 MeV for J/ipK%, \AE\ < 30 MeV 
for iP{2S)K% and \AE\ < 25 MeV for Xc\K%) as well as 
the p* B distribution for J/ipK® candidates are shown in 
Fig. [TJ We require 5.27 GeV/c 2 < M hc < 5.29 GeV/c 2 
for fcp modes with a K° s and 0.20 GeV/c < p* B < 
0.45 GeV/c for J/^K° L for the fit to the CP violation 
parameters. For the candidates passing all the criteria 
mentioned above, the signal yield and purity are esti- 
mated for each CP-eigenstate and listed in Table [I] 

We determine Sf and Af for each mode by perform- 
ing an unbinned maximum-likelihood fit to the observed 
At distribution. The probability density function (PDF) 
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FIG. 1: (color online), (a) Mhc distribution within the AE 
signal region for B° ->■ J/tpK% (black), tp{2S)K° s (blue), and 
XciKg (magenta); the superimposed curve (red) shows the 
combined fit result for all these modes, (b) p* B distribution of 
B° — > J/rpKl candidates with the results of the fit separately 
indicated as signal (open histogram), background with a real 
J/i/j and real K^'s (yellow), with a real J/ip and a fake K% 
candidate (green), and with a fake J/tp (blue). 



TABLE I: CP eigenvalue (£/)> signal yield (-/V s i g ) and purity 
for each B° — > fcp mode. 



Decay mode 



i>{2S){t+l-)K° s 
iP(2S)(J/^+ty-)K° s 

J/i>Kl 



Purity (%) 



-1 12649±114 

-1 904± 31 

-1 1067± 33 

-1 940± 33 

+1 10040±154 



97 
92 
90 
86 
63 



for the signal distribution, P s ; g (At; Sf,Af,q,wi, Awi), is 
given by Eq. (TO. fixing t b o and Am^ at their world 
average values |16j and including modifications to take 
the effect of incorrect flavor assignment (parameterized 
by wi and Awi) into account. The distribution is con- 
volved with the proper-time interval resolution function, 
Psig(Ai), formed by convolving four components: the de- 
tector resolutions for zcp and ztag, the shift of the zt ag 
vertex position due to secondary tracks from charmed 
particle decays, and the kinematic approximation that 
the B mesons are at rest in the CM frame [17[. Because 
we now use h to characterize the vertex goodness-of-fit, 
each of these resolution function components in Ref. [l~7| 
is reformulated as a function of h and a z . 

Using the M^ c sideband events, the background PDF, 
Pbkg(Ai), for each of the CP-odd modes is modeled as 
a sum of exponential and prompt components, and is 
convolved with Rhkg(At) expressed as a double Gaussian 
function. In the J /ipK^ mode, there are CP violating 
modes among the B — > J/tjjX backgrounds, which are 
included in the background PDF. The At PDFs for the 
remaining B — > J/ipX and other combinatorial back- 
grounds are estimated from the corresponding large MC 
sample and M^+ e ~ sideband events, respectively. The 
construction of these PDFs follows the same procedure 
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as in our previous analyses @, [HJ . 

We determine the following likelihood for the i-th 
event: 

Pi = (1- / [PkiA^RkiAU - At')]d(At') 



-f ol P ol (AU), 



(2) 



where the index k labels each signal or background com- 
ponent. The fraction fk depends on the r region and 
is calculated on an event-by-event basis as a function 
of AE and Mb c for the CP-odd modes and p* B for the 
CP-even mode. The term P Q \(At) is a broad Gaussian 
function that represents an outlier component / i, which 
has a fractional normalization of order 0.5% [17[ . The 
only free parameters in the fits are Sf and Af, which 
are determined by maximizing the likelihood function 
L = J\ i Pi(Ati;Sf,Af). This likelihood is maximized 
for each fcp mode individually, as well as for all modes 
combined taking into account their CP-eigenstate val- 
ues; the results are shown in Table [TT1 Figure [5] shows the 
At distributions and asymmetries for good tag quality 
(r > 0.5) events. We define the background-subtracted 
asymmetry in each At bin by (N + — N-)/(N+ + AC), 
where AC (AC) is the signal yield with q = +!(—!). 




-6 -4 -2 2 4 6 

At (ps) 



-6 -4 -2 2 4 6 
At (ps) 



FIG. 2: (color online) The background-subtracted At distri- 
bution (top) for q = +1 (red) and q = — 1 (blue) events and 
asymmetry (bottom) for good tag quality (r > 0.5) events 
for all CP-odd modes combined (left) and the CP-even mode 
(right). 



TABLE II: CP violation parameters for each B° —> fcp mode 
and from the simultaneous fit for all modes together. The first 
and second errors are statistical and systematic uncertainties, 
respectively. 



TABLE III: Systematic errors in Sf and Af in each fcp mode 
and for the sum of all modes. 



Decay mode sin20i = — 5/5/ 



i>{2S)K% 
J/i>K L 



Af 



+0.670 ± 0.029 ±0.013 
+0.738 ± 0.079 ± 0.036 
±0.640 ±0.117 ±0.040 
+0.642 + 0.047 + 0.021 



-0.015 + 0.021 
+0.104 + 0.055 
-0.017 + 0.083 
+0.019 + 0.026 



+0.045 

0.023 
+0.047 
0.027 
+0.046 
0.026 
+0.017 
0.041 



All modes +0.667 ± 0.023 ± 0.012 ±0.006 ± 0.016 ± 0.012 



Uncertainties originating from the vertex reconstruc- 
tion algorithm are a significant part of the systematic 
error for both sin20i and Af. These uncertainties are 
reduced by almost a factor of two compared to the previ- 
ous analysis [7] by using h for the vertex-reconstruction 
goodness-of-fit parameter, as described above. In partic- 
ular, the effect of the vertex quality cut is estimated by 
changing the requirement to either h < 25 or /i < 100; the 
systematic error due to the IP constraint in the vertex re- 
construction is estimated by varying the IP profile size in 
the plane perpendicular to the z-axis; the effect of the cri- 
terion for the selection of tracks used in the /t ag vertex is 
estimated by changing the requirement on the distance of 
closest approach with respect to the reconstructed vertex 
by ±100 /im from the nominal maximum value of 500 /im. 
Systematic errors due to imperfect SVD alignment are es- 
timated from MC samples that have artificial misalign- 
ment effects. Small biases in the Az measurement are 
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±0.003 


±0.003 


±0.003 


Possible 


s f 


±0.004 


±0.004 


±0.004 


±0.004 


±0.004 


fit bias 


Af 


±0.005 


±0.005 


±0.005 


±0.005 


±0.005 


Signal 


Sf 


±0.004 


±0.016 


< 0.001 


±0.016 


±0.004 


fraction 


Af 


±0.002 


±0.006 


< 0.001 


±0.006 


±0.002 


Background 


Sf 


< 0.001 


±0.002 


±0.030 


±0.002 


±0.001 


At PDFs 


Af 


< 0.001 


< 0.001 


±0.014 


< 0.001 


< 0.001 


Physics 


s f 


±0.001 


±0.001 


±0.001 


±0.001 


±0.001 


parameters 


Af 


< 0.001 


< 0.001 


±0.001 


< 0.001 


< 0.001 


Total 


s f 


±0.013 


±0.036 


±0.040 


±0.021 


±0.012 




Af 


+0.045 
-0.023 


+0.047 
-0.027 


+0.046 
-0.026 


+0.017 
-0.041 


±0.012 



observed in e + e~ — > /C/x - and other control samples: to 
account for these, a special correction function is applied 
and the variation with respect to the nominal results is 
included as a systematic error. We also vary the \At\ 
range by ±30 ps to estimate the systematic uncertainty 
due to the \At\ fit range. The vertex resolution function 
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is another major source of sin20i and Af uncertainty. 
This effect is estimated by varying each resolution func- 
tion parameter obtained from data (MC) by ±ler (±2<r) 
and repeating the fit to add each variation in quadrature. 
The uncertainty in the estimated errors of the parameters 
of reconstructed charged tracks is also taken into account. 
The largest contribution to the systematic uncertainty in 
Af is the effect of the tag- side interference (TSI), which 
is described in detail in [18[ . Since the effect of TSI has 
opposite sign for different CP-eigenstates, there is a par- 
tial cancellation in the combined result. Hence the com- 
bined TSI systematic is smaller than the systematic in 
each individual mode. Systematic errors due to uncer- 
tainties in the wrong-tag fractions are studied by varying 
the wrong-tag fraction individually in each r region. A 
possible fit bias is examined by fitting a large number of 
MC events. Other contributions come from uncertainties 
in the signal fractions, the background At distribution, 
tbo and Amd- Each contribution is summarized in Ta- 
ble IIIII We add them in quadrature to obtain the total 
systematic uncertainty. 

In summary, we present the final sin2(/>i measure- 
ment using the entire Belle T(45) data sample con- 



taining 772 x 10 6 BB pairs. We have reconstructed 
b — > ccs induced B meson decays in three CP-odd 
modes {J/^K%, ^{2S)K%, and XdK° s ) and one CP- 
even mode (J/ipK^). The fit, using common CP- 
sensitive parameters for all four modes, yields the val- 
ues sin 2<fo = 0.667 ±0.023(stat)±0.012(syst) and A f = 
0.006 ± 0.016(stat) ± 0.012(syst). The results are consis- 
tent with previous measurements @, Q ■ These are the 
most precise determination of these parameters and so- 
lidify the SM reference value used to test for evidence of 
new physics beyond the SM. 
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